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Abstract: Results for next-to-leading order QCD corrections to the pp(pp) — > tt — > 
W + W~bb — > e + v e ^~v^bb +X processes with complete off-shell effects are presented 
for the first time. Double-, single- and non-resonant top contributions of the order 
0{a z s a i ) are consistently taken into account, which requires the introduction of a 



complex-mass scheme for unstable top quarks. Moreover, the intermediate W bosons 
are treated off-shell. Comparison to the narrow width approximation for top quarks, 
where non-factorizable corrections are not accounted for is performed. Besides the 
total cross section and its scale dependence, several differential distributions at the 
TeVatron run II and the LHC are given. In case of the TeVatron the forward- 
backward asymmetry of the top is recalculated afresh. With inclusive selection cuts, 
the forward-backward asymmetry amounts to A l FB = 0.051 ± 0.0013. Furthermore, 
the corrections with respect to leading order are positive and of the order 2.3% for 
the TeVatron and 47% for the LHC. A study of the scale dependence of our NLO 
predictions indicates that the residual theoretical uncertainty due to higher order 
corrections is 8% for the TeVatron and 9% for the LHC. 
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1. Introduction 

The ti production process is a copious source of W^-pairs and, hence, of isolated 
leptons at the TeVatron and the LHC. In consequence it is intensely studied as a 
signal at these colliders. In view of the large production rate, precise and direct 
measurements are possible, which require a detailed theoretical understanding. In 
addition, it constitutes an important background for many new particle searches. 
Examples include the leptonic signals for cascade decays of supersymmetric particles 
or searches for H — » W + W and H — > t + t~ decays. 

Even though, the first results for next-to-leading order (NLO) QCD corrections 
to heavy quark production were presented more than twenty years ago [1-4], recent 
progress in NLO [5-9] and next-to-next-to leading order (NNLO) [10-15] calculations, 
as well as next-to-next-to-leading-log resummations (NNLL) [16-22] for inclusive ti 
hadroproduction is truly astonishing. 

The list for the more exclusive channels is just as impressive: NLO QCD correc- 
tions have been calculated for the tiH signal [23-28] , where the Higgs boson has been 
treated as a stable particle. Most recently the factorizable QCD corrections to this 
process have been presented [29], where higher order corrections to both production 
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and decay of the Higgs boson to a bb pair have been calculated with the latter mod- 
eled by the Higgs propagator with a fixed width. Moreover, NLO QCD corrections to 
a variety of 2 — > 3 backgrounds processes ttj [30-32] , tiZ [33] and ttj [34] have been 
obtained. Most recently, NLO QCD corrections to 2 — > 4 backgrounds tibb [35-38] 
and tijj [39] have also been evaluated. 

Usually, ti production is restricted to on-shell states and decays if available are 
treated in the narrow- width approximation (NWA), which effectively decouples top 
production and decay. The NWA allows to neglect non-resonant as well as non- 
factorizable amplitude contributions, thus leading to significant simplifications for 
calculations of higher order corrections. Whenever resonant top production domi- 
nates, as it does for very inclusive cuts, this approximation is of course well motivated. 
In some cases calculations have been further simplified by also treating the decaying 
W bosons as on-shell particles. 

Naturally, the accuracy of these approximations needs to be tested, which re- 
quires a full calculation of off-shell effects. One thus needs a calculation which in- 
cludes both resonant and non-resonant contributions, using finite width top-quark 
propagators, which correctly includes interference effects between the various contri- 
butions. The purpose of this paper is to present such a complete calculation for ti 
production at NLO QCD level. In addition to merging resonant and non-resonant 
effects for the top quarks, we also include finite width effects for the W bosons, i.e. 
we consider NLO QCD corrections to the general e + v e ^L~v^bb final state. 

In addition, all selection strategies based on next-to-leading order simulations, 
which have been devised for the efficient suppression of ti background, are at present 
optimized against top production in the NWA. Within our approach, presented in the 
form of a flexible Monte Carlo program which allows to study NLO QCD corrections 
to cross sections and kinematic distributions with arbitrary cuts on particles in the 
final state and with full spin correlations, it is possible to reexamine the quality of 
the chosen selection with improved accuracy. 

The paper is organized as follows. In Section |2| we briefly describe the calculation 
of the NLO corrections. Numerical results for the integrated and differential cross 
sections are presented in Section |3] both for the TeVatron and the LHC. Finally, we 
conclude in Section f|. 

2. Theoretical framework 
2.1 Born level 

At Born level the partonic reactions are 

gg e + v e ^v il bb 
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Figure 1: Representative Feynman diagrams contributing to the leading order process 
gg — > e ir v e yTv ll bb at 0(a 2 s cfi), with different off-shell intermediate states: double-, single-, 
and non-resonant top quark contributions. 

qq -)■ e + u e ^u fl bb (2.1) 

where q stands for up- or down-type quarks. The 0{a 2 s a 4 ) contributions to the 
e + v e ^v^bb process can be subdivided into three classes, namely diagrams containing 
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two top quark propagators that can become resonant, diagrams containing only 
one top quark resonance and finally diagrams without any top quark resonance. 
Regarding the W ± resonances one can distinguish only two subclasses, double- and 
single-resonant gauge boson contributions. A few examples of Feynman diagrams 
contributing to the leading order gg — > e + u e fi~u fJi bb subprocess are presented in 
Figure [I]. 

Since the produced top quarks are unstable particles, the inclusion of the decays 
is performed in the complex mass scheme, which for LO is described in Ref. [40,41]. 
It fully respects gauge invariance and is straightforward to apply. In the amplitude 
(at LO and NLO) we simply perform the substitution 

-m t -Me) -1 (rf-fit + i4 = nit - im t r t . (2.2) 

Since we are interested in NLO QCD corrections, gauge bosons are treated within 
the fixed width scheme. Our LO results have been generated with the Helac- 
Dipoles [42] package and cross checked with Helac-Phegas [43,44], a generator 
for all parton level processes in the Standard Model, which has, on its own, already 
been extensively used and tested in phenomenological studies see e.g. [45-49]. The 
integration over the fractions x\ and of the initial partons is optimized with 
the help of Parni [50]. The phase space integration is executed with the help of 
Kaleu [51] and cross checked with Phegas [52], both general purpose multi-channel 
phase space generators. 

Furthermore, results have been checked against another program that computes 
the ti production cross section with top decays, namely Mcfm [53]. A perfect 
agreement has been found with our results, both for the TeVatron and the LHC, 
once top quarks and W gauge bosons have been put on shell in the Helac-Dipoles 
package. We additionally reproduced results presented in Ref. [8] again assuming 
that both tops and W's are on shell. 

2.2 The virtual corrections 

The virtual corrections consist of the 1-loop corrections to the LO reactions. One 
can classify the corrections into self-energy, vertex, box-type, pentagon-type and 
hexagon-type corrections. Typical examples of the virtual graphs are shown in Fig- 
ure [|. In evaluating the virtual corrections, the Helac-ILoop [54] approach is 
used. It is based on the Helac-Phegas program to calculate all tree-order like 
ingredients and the OPP [55] reduction method. The cut-constructible part of the 
virtual amplitudes is computed using the CutTools [56] code. The rational term 
Ri of the amplitude is computed by the CutTools code as well, whereas the R2 
term, by the use of extra Feynman rules as described in [56,57]. Numerical results 
are obtained using the same methods as described in [37]. 
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Figure 2: Representative Feynman diagrams contributing to the virtual corrections to the 
partonic subprocess gg — > e + i' e ^i~i>^bb at 0(a^a 4 ). 

As explained before, the process under consideration requires a special treatment 
of unstable top quarks, which is achieved within the complex- mass scheme [40]. At 
the one-loop level the appearance of a non-zero top-quark width in the propagator 
requires the evaluation of scalar integrals with complex masses, for which the program 
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OneLOop [54, 58] is used. We also need mass renormalization for the top quark, 
which, for consistency, is done by using a complex mass in the well known on- 
shell mass counterterm. The preservation of gauge symmetries (Ward Identities) 
[40,59-61] by this approach has been explicitly checked up to the one-loop level. 

Although finite width effects have been studied routinely at tree order, the same 
is not true for calculations at the one loop level. A novel aspect of the introduc- 
tion of a non-zero width is the effect on the infrared structure of the scattering 
amplitudes. Working in dimensional regularization, soft and collinear singularities 
arise. When massive particles acquire a complex mass, the soft 1 /e-singularities 
due to the exchange of gluons, are replaced by factors proportional to log(T t /m t ), 
that become singular in the limit T t — > 0. We have explicitly checked that in- 
cluding all contributions, factorizable and non-factorizable, the usual cancellation of 
infrared 1/e 2 and 1/e poles between virtual and real corrections, the latter repre- 
sented by the X(e)-operator, takes place. This means that a partial cancellation of 
\og(T t /m t ) terms happens within the virtual corrections alone. Nevertheless loga- 
rithmic enhancements remain in the finite part of the virtual corrections and have to 
be cancelled by corresponding terms from the real corrections, since they represent 
the same soft singularities, dimensionally regularized in the case of on-shell particles. 



2.3 The real emission 

The generic processes for the real corrections are given by 

gg e + v e ix~V lx bbg 
qg -» e + v e n-v^bbq 

gq e + u e fi'u fl bbq (2.3) 
qq e + u e fi~u fl bbg 

(where again q stands for up- or down-type quarks) and include all possible contri- 
butions of the order of 0(a^a 4 ). The complex mass scheme for unstable top quarks 
has been implemented in complete analogy to the LO case. 

We employ the dipole subtraction formalism [62] to extract the soft and collinear 
infrared singularities and to combine them with the virtual corrections. Specifically, 
the formulation [63] for massive quarks has been used with the extension to arbitrary 
helicity eigenstates of the external partons [42], as implemented in Helac-Dipoles. 
In the case at hand, the number of dipoles is as follows: 27 for the process gg — > 
e^VziTv^bbg and 15 for processes qq — > e^u^u^bbg, qg — > e^u^u^bbq and gq — > 
e + u e fi~i? fl bbq. Let us stress at this point, that, similarly to most authors, we do 
not use finite dipoles regularizing the quasi-collinear divergence induced by both 
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top quarks moving in the same direction, even though they are implemented in the 
software. Due to the large top quark mass, they do not improve numerical stability. 

Besides the cancellation of divergences, which we have mentioned in the previous 
section, we have also explored the independence of the results on the unphysical cutoff 
in the dipole subtraction phase space (see [42] and references therein for details) to 
further check our calculation. 

2.4 Phase space generation 

In LO calculations, the jet definition consists of a set of phase space cuts not allowing 
any parton to become arbitrarily soft, and no pair of partons to become arbitrar- 
ily collinear. This changes for the real-radiation contribution in NLO calculations, 
for which single partons are allowed to become arbitrarily soft and single pairs of 
partons are allowed to become arbitrarily collinear. This means that phase space 
generators like PHEGAS [52] and KALEU [51], which construct momentum configu- 
rations from kinematical variables generated following a priori defined probability 
densities, cannot be directly applied in their LO set-up, since these densities antici- 
pate the singular behavior of the squared amplitudes, and are typically not defined 
in the soft and collinear limits. Furthermore, the subtraction terms in the dipole- 
subtraction scheme, used to eliminate the singularities in the real-radiation phase 
space integral, do not exactly follow the same peak structure as the tree-level n + 1- 
particle matrix element squared, whereas Phegas and Kaleu are designed only to 
efficiently deal with the latter. We chose to deal with this situation via a multi- 
channel approach [64], in which a separate channel is associated with each term in 
the real-subtracted integral, i.e., with the tree-level n + 1-particle matrix element 
squared as well as each dipole term. 

The channel for the n + 1-particle matrix element squared generates momenta 
using an instance of Kaleu anticipating the peak structure of this integrand. The 
phase space defined by promoting the LO cuts to n + 1 partons is filled in the usual 
LO approach. The soft and collinear regions "below the cuts" are filled by replacing 
the densities for the invariants by densities that are integrable in these regions. 

All dipole channels also carry their own instances of Kaleu, but each of these 
generates n-momentum configurations anticipating the peak structure of the n- 
particle matrix element squared of the underlying process of the dipole term. Such a 
n-momentum configuration is then turned into an n + 1-momentum configuration by 
essentially applying the inverse of the phase space mapping performed in the calcula- 
tion of the dipole contribution itself. This generation of an extra momentum follows 
exactly the formulas for the parton showers based on the dipole formalism presented 
in [65] and [66] . The azimuthal angle needed for the construction of the extra momen- 
tum is generated with a flat distribution, and the other two variables, traditionally 
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denoted (yij,k,Zi) for final-final, (xij ta ,Zi) for final-initial, (xij j0/ ,Ui) for initial-final, 
and (xi ya b, Vi) for initial-initial dipoles, are generated following self-adaptive densi- 
ties. This happens "on the fly" during the Monte Carlo integration, following the 
approach presented in [50]. Finally, each instance of Kaleu carries a mult i- channel 
weight in the "highest level" multi-channel density which is optimized during the 
Monte Carlo integration, and each instance performs its own internal multi-channel 
optimization, as described in [51]. 

We have performed a few tests to check the performance of this new approach in 
case of the qq — > e + u e ji^u^bbg subprocess. More precisely, we have made a compar- 
ison between three options, namely Kaleu with dipole channels, Kaleu without 
dipole channels and Phegas, which does not have dipole channels. Since the compu- 
tational cost comes mainly from the accepted events, comparisons are made at equal 
numbers of accepted events. Our findings can be summarized as follows. For the 
dipole phase space cut-off parameter a max = 1, when all dipoles are calculated for 
each phase space point, Phegas and Kaleu without dipoles channels are compa- 
rable in terms of errors. KALEU with dipole channels, however, gives an error which 
is 5 times smaller. Realize that this implies a reduction in the number of events 
by a factor of 25 to reach the same error. For a max = 0.01, when much less dipole 
subtraction terms are needed per event, the improvement is not so dramatic, and 
the introduction of the dipole channels reduces the error by a factor 3 compared to 
Kaleu without dipole channels, and a factor 2 compared to Phegas, implying a 
reduction in necessary events by a factor of 9 and 4 respectively. 

We conclude that the dipole channels structurally improve the convergence of 
the phase space integrals for the real-subtracted contribution. It is, however, difficult 
to express the improvement quantitatively because it depends on the process and the 
value of parameters like a max . 

3. Numerical Results 
3.1 Setup 

We consider the process pp(pp) — > ti + X — > W + W~bb + X — > e + u e fi^i , ^bb + X both 
at the TeVatron run II and the LHC i.e. at a center-of-mass energy of s/s = 1.96 TeV 
and \fs = 7 TeV correspondingly. For the LHC case we additionally calculate the 
integrated cross section at a center-of-mass energy y/s = 10 TeV. We only simulate 
decays of the weak bosons to different lepton generations to avoid virtual photon 
singularities stemming from quasi-collinear 7* — > £ + £~ decays. These interference 
effects are at the per-mille level for inclusive cuts, as checked by an explicit leading 
order calculation. The complete ifi^ cross section (with = e, fi) can be obtained 
by multiplying the result with a lepton-flavor factor of 4. We keep the Cabibbo- 
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Kobayashi-Maskawa mixing matrix diagonal. The unstable (anti)top quark is treated 
within the (gauge invariant) complex-mass scheme, as explained in the previous 
section. The Standard Model parameters are given the following values within the 
scheme [67]: 

m w = 80.398 GeV, Y w = 2.141 GeV, 
m z = 91.1876 GeV, T z = 2.4952 GeV, 
= 1.16639 x 1(T 5 GeV" 2 , 

sin 2 6w = 1 — r mw/ m 'z- (3-1) 
The electromagnetic coupling is derived from the Fermi constant G M according to 

^ = V^Gf.m^ sin 2 W ^ 

7T 

For the top quark mass we take m t = 172.6 GeV and all other QCD partons including 
6 quarks as well as leptons are treated as massless. The contribution from the Higgs 
boson can be neglected since for inclusive cuts it is below 1%. In our case, however, 
the 6-quarks are massless and the Higgs contribution simply vanishes. The top quark 
width calculated from [68,69] is = 1.48 GeV at LO and T? LO = 1.35 GeV at 
NLO where a s = a s (m t ) = 0.107639510785815. Mass renormalization is performed 
in the on-shell scheme. All final-state 6 quarks and gluons with pseudorapidity \r)\ < 5 
are recombined into jets with separation yj A<p 2 + Ay 2 > D = 0.4 in the rapidity- 
azimuthal angle plane via the following IR-safe algorithmes: the kx algorithm [70-72] , 
the anti-kr algorithm [73] and the inclusive Cambridge/ Aachen algorithm (C/A) [74]. 
The distance measure dij for these algorithms is defined as 

AR 2 

= min {p%(i),p%(j)) ~JyT 
d iB =PT(i) , 



where AR^ = w A<pf - + Ay 2 - and the parameter p is equal to 1 for the kx algorithm, 

for C/A and —1 for anti-kx algorithm. Moreover, we impose the following additional 
cuts on the transverse momenta and the rapidity of two recombined 6-jets: pr b > 20 
GeV, \y b \ < 4.5 where 



Basic selection is applied to (anti)top decay products to ensure that the leptons 
are observed inside the detector and are well separated from each other: px e > 
20 GeV, \rj e \ < 2.5, AR je > 0.4, where j = 6,6, and p Tmias > 30 GeV. In the 
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following we consistently use the CTEQ6 set of parton distribution functions (PDFs) 
[75,76]. More precisely, we take CTEQ6L1 PDFs with a 1-loop running a s in LO and 
CTEQ6M PDFs with a 2-loop running a s in NLO. The contribution from b quarks 
in the initial state is neglected, since at LO for inclusive cuts this contribution is 
suppressed to the per-mille level. The number of active flavors is Np = 5, and 
the respective QCD parameters are = 165 MeV and A^' IS = 226 MeV. In the 
renormalization of the strong coupling constant, the top-quark loop in the gluon 
self-energy is subtracted at zero momentum. In this scheme the running of a s is 
generated solely by the contributions of the light-quark and gluon loops. By default, 
we set the renormalization and factorization scales, fiR and /xp, to the common 
value fj, = m t . For inclusive cuts, where the contribution from the double resonance 
Feynman diagrams dominates, the top mass is a valid scale. 

3.2 Results for the TeVatron run II 

We begin our presentation of the final results of our analysis with a discussion of the 
total cross section at the central value of the scale, /ij? = fip = m t at the TeVatron 
run II. The respective numbers are presented in Table [l] for the two choices of the 
dipole phase space cutoff parameter a max (see e.g. [42] for more details) and for three 
different jet algorithms. At the central scale value, the full cross section receives small 
NLO correction of the order of 2.3%. 

Subsequently, we turn our attention to the scale dependence for the total cross 
section at LO and NLO. The left panel of Figure |3] shows the dependence of the 
integrated LO cross section on the renormalization and factorization scales where 
A* — /Ar — fJ-F — £, m t- The variation range is taken from \i = m t /8 to \i = 8m t . 
The dependence is large, illustrating the well known fact that the LO prediction 
can only provide a rough estimate. At the TeVatron with our cut selection the qq 
channel dominates the total pp cross section by about 95% followed by the gg channel 
with about 5%. In the right panel the scale dependence of the NLO cross section is 
shown together with the LO one. As expected, we observe a reduction of the scale 
uncertainty while going from LO to NLO. Varying the scale down and up by a factor 
2 changes the cross section by +40% and —26% in the LO case, while in the NLO 
case we have obtained a variation of the order —8% and —4%. Let us mention here 
that while calculating the scale dependence for the NLO cross section we kept Yf 
fixed independently of the scale choice. The error introduced by this treatment is 
however of higher order, and particularly for two scales /i = m t /2 and \x = 2mt 
amounts to ±1.5% respectively. 

In the following we would like to estimate the size of the non-factorizable correc- 
tions for our inclusive setup. To achieve this the full result has been compared with 
the result in the NWA. The latter has been obtained by rescaling the coupling of 
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Algorithm 




a-£o =1 [fb] 


<lo =001 [fb] 


antirhr 


34.922 ± 0.014 


35.705 ± 0.047 


35.697 ± 0.049 


kr 


34.922 ± 0.014 


35.727 ± 0.047 


35.723 ± 0.049 


C/A 


34.922 ± 0.014 


35.724 ± 0.047 


35.746 ± 0.050 



Table 1: Integrated cross section at LO and NLO for pp — > e + u e fj 1 ~P fl bb + X production 
at the TeVatron run II with ^/s = 1.96 TeV, for three different jet algorithms, the anti-kx, 
kx and the Cambridge/ Aachen jet algorithm. The two NLO results refer to different values 
of the dipole phase space cutoff a max - The scale choice is fj,R = (j,f = rut- 
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Figure 3: Scale dependence of the LO cross section with the individual contributions of 
the partonic channels (left panel) and scale dependence of the LO and NLO cross sections 
(right panel) for the pp —> e + v e ^~v^bb + X process at the TeVatron run II with y/s = 1.96 
Te V, where renormalization and factorization scales are set to the common value \i = fj,R = 



the top quark to the W boson and the b quark by several large factors to mimic the 
limit T t — y when the scattering cross section factorizes into on-shell production and 
decay. Our findings are depicted in Figure f| where the dependence of the total NLO 
cross section together with its individual contributions, real emission part and LO 
plus virtual corrections, are shown. The behavior is compatible with a logarithmic 
dependence on T t , which cancels between real and virtual corrections. For inclu- 
sive production, advancing from NWA to the full result changes the cross section no 
more than +1% which is consistent with the uncertainty of the NWA i.e. of order 
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Figure 4: Dependence of the NLO cross section, cjt, (red solid line) and the individual 
contributions, the real emission part, cxr, (green dashed line) and the LO plus virtual part, 
ay, (blue dotted line), on the rescaling parameter ( defined as T resca [ ec i = (T t for the 
pp — > e + v e ^~V[jbb + X process at the TeVatron run II with y/s = 1.96 TeV. Act is defined 
as follows: Aa^C) = (^(C) - o"i(C = l))/cr T (C = 1) with i = V,R,T. 

0(T t /m t ). 

Comparing our NLO integrated cross section with the value cxnlo = 36.47 fb 
presented in Ref . [8] , we observe a 2% discrepancy, which can easily be explained by 
two effects. First of all, in [8] NLO QCD corrections have been calculated employing 
an on-shell approximation for the top quarks and the W bosons. The former approx- 
imation can introduce a difference of the order of 0(T t /mt) ~ 1% while the latter 
of the order of 0(Y W /m w ) ~ 3%. As a second effect, there are small differences 
between individual setups, in e.g. the value of T t , m t , PT miS3 and ARj£. 

We have also compared our results with those generated with Mcfm. We have 
been able to use the same cuts and input parameters, but there is an essential 
difference as far as the construction of the cross section is concerned. Indeed, Mcfm 
includes corrections to the production of on-shell top quarks only, whereas decays are 
included at leading order. Moreover, W bosons are also treated in the narrow width 
approximation. In the end, Mcfm gives the following results ctlo = (36.494 ± 0.050) 
fb and unlo = (39.622 ± 0.065) fb, which are different from ours by 4.5% at LO and 
by 11% at NLO. Although we have not quantified the impact of different approaches 
used, related to the top quark and the W-boson finite width, as well as the NLO 
corrections to the decay of the top quarks, the overall comparison seems reasonable 
and compatible with estimates based on the order of magnitude for these effects, 
0(T t /m t ), 0(Tw l m w) an d 0(a s ). A more detailed study would be necessary in 
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Figure 5: Differential cross section distributions as a function of rapidity yt of the top 
(red solid curve) and anti-top quarks (blue dotted curve), rapidity yi, of the b-jet (red solid 
curve) and anti-b-jet (blue dotted curve) and rapidity yi of the positron (red solid curve) and 
muon (blue dotted curve) at next-to-leading order for the pp — > e + v e fj,~v , fJj bb + X process 
at the TeVatron run II. The green dashed curves correspond to the leading order results. 

order to establish the relevance of these differences for the experimental analysis, 
which goes beyond the purpose of the present publication. 

In a next step we recalculate the top quark forward-backward asymmetry for the 
TeVatron from the top rapidity distribution. We show our results for the LO and 
NLO inclusive calculations. At LO, tt production is totally charge-conjugation sym- 
metric for both production mechanisms (quark and gluon fusion). As a consequence, 
the angular distributions of the t and i are symmetric with respect to the beam 
axis for pp collisions. However, at higher orders in a s , this is not longer true. Not 
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Figure 6: Differential cross section distributions as a function of the invariant mass m t i 
of the top-anti-top pair, rapidity y t i of the top-anti-top pair, averaged transverse momentum 
PT t of the top and anti-top and averaged rapidity yt of the top and anti-top for the pp —> 
e + v e ii~v^bh + X process at the TeVatron run II. The blue dashed curve corresponds to 
the leading order, whereas the red solid one to the next-to-leading order result. The lower 
panels display the differential K factor. 

all processes involving additional partons are symmetric under charge conjugation 
with respect to the incoming parton and anti-parton beams. As was pointed out in 
Ref. [77-79] the process gg —> tig is, but the processes qq — >■ tig and qg — > tiq are not. 
Processes involving initial state valence quarks will therefore exhibit a charge asym- 
metry. This is caused by interference between initial and final state gluon emission 
on the one side and by interference between color singlet 4-point virtual corrections 
and the Born term for the qq process [1,3] on the other. Because it production at 
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Figure 7: Differential cross section distributions as a function of the averaged transverse 
momentum px b of the b-jet and anti-b-jet, averaged rapidity yb of the b-jet and anti-b- 
jet and AR b i separation for the pp — > e + v e p^U^bb + X process at the TeVatron run II. 
The blue dashed curve corresponds to the leading order, whereas the red solid one to the 
next-to -leading order result. The lower panels display the differential K factor. 



the TeVatron is dominated at the 95% level by qq annihilation, as was mentioned 
earlier in the paper, we can expect the qq subprocess asymmetry to be visible in the 
total sample. The integrated charge asymmetry is defined through 

, 4>o^(y)-4 >0 iW 

i>oW + 4>o^/)' 

where y t (yi) is the rapidity of the top (anti-top) quark in the laboratory frame and 
N t (y) = da t t/dy t , Nt(y) = da t i/dyi. Due to the CP invariance of QCD the rapid- 
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Figure 8: Differential cross section distributions as a function of the averaged transverse 
momentum pr w of the bosons and averaged rapidity yw of the bosons for the 
pp —> e + v e fi~v i l j j bb -\- X process at the TeVatron run II. The blue dashed curve corresponds 
to the leading order, whereas the red solid one to the next-to-leading order result. The lower 
panels display the differential K factor. 

ity distributions of top and anti-top are mirror images of each other, i.e. iVj(y) = 
N t (—y), and integrated charge asymmetry is equal to the integrated forward-backward 
asymmetry of the top quark defined as 

Moreover, A FB = —A t FB . 

As can be seen in the upper-left part of the Figure |5] the LO ti inclusive cross 
section is symmetric around yt — (green dashed curve). The NLO inclusive result 
for the top/anti-top quark is, on the other hand, shifted to larger y t for the top 
quark (solid red curve) and smaller y t for the anti-top quark (dotted blue curve). 
This corresponds to a positive integrated forward-backward asymmetry of the order 
of 

A FB = 0.051 ±0.0013, (3.5) 

which tells us that top quarks are preferentially emitted in the direction of the in- 
coming protons. 

Next-to-leading order contributions to the forward-backward asymmetry have 
already been calculated in the on-shell ti production [80] and amount to A FB = 
0.051 ±0.006. The CDF measurement based on 5.3 fb^ 1 integrated luminosity in the 



-16- 





Figure 9: Differential cross section distributions as a function of the averaged transverse 
momentum pT e of the charged leptons, averaged rapidity yi of the charged leptons, PT misa 
and ARu for the pp — > e Jr v e \x~v iJ bb + X process at the TeVatron run II. The blue dashed 
curve corresponds to the leading order, whereas the red solid one to the next-to-leading 
order result. The lower panels display the differential K factor. 



semi-leptonic channel yields A FB = 0.150 ± 0.050 s * a *- ± 0.02A syst [81], while the D0 
measurement of this asymmetry yields A FB = 0.08 ±0.04 sto *' ±0.01 s ^*' based on 4.3 
fb _1 integrated luminosity [82]. The uncertainties of these results are still very large 
and statistically dominated. 

In the same manner we can calculate the integrated forward-backward asymme- 
try for the top decay products, namely the 6-jet and the positively charged lepton. 
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Figure 10: Differential cross section distribution as a function of the total transverse 
energy, Ht, for the pp — > e + u e fi~u f j i bb +X process at the TeVatron run II. The blue dashed 
curve corresponds to the leading order, whereas the red solid one to the next-to-leading order 
result. The lower panels display the differential K factor. 

Our results can be summarized as follows: 

A b FB = 0.033 ± 0.0013 , A l p B = 0.034 ± 0.0013 , (3.6) 

where 

Ab ¥+ _ 4 + >o^)-4 + <o^) ^ 

FB 4>o^(^) + 4<oW FB I yi+ >o^ + (y) + I yi+<0 Ne + (y) [6J) 

and ui and y^ are the rapidity of the charged lepton and the 6-jet respectively and 
Ni+{y) = d<j t t/dyi+, Nb(y) = doiijdy^. In case of A e FB we agree with Ref. [7] where 
A e FB = 0.033 has been quoted. The integrated forward-backward asymmetries of the 
charged lepton and the 6-jet have the same sign as A FB but are smaller in magnitude. 
Let us stress at this point, that the 6-jet integrated forward-backward asymmetry is 
a rather theoretical observable even though it can in principle be measured once the 
6-jet is distinguished experimentally from the anti-6-jet through e.g. the charge of 
the associated lepton flying in the same direction. However, it is extremely difficult 
to determine the charge of the 6-jet and this measurement will heavily depend on 
the 6-jet tagging efficiency The 6-jet and charged lepton differential distributions in 
rapidity are also presented in Figure || 

While the size of the corrections to the total cross section is certainly interesting, 
it is crucial to study the corrections to distributions. In the following, the NLO QCD 
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corrections to the differential distributions for the dileptonic channel with full off-shell 
effects are presented. 

In Figure || we start with the most important observable, namely, the differential 
distribution of the tt invariant mass, m t i . Figure |] depicts also the rapidity, y t t, of 
the top-anti-top system as well as the averaged transverse momentum, pT t , and the 
averaged rapidity y t of the top and anti-top. The blue dashed curve corresponds 
to the leading order, whereas the red solid one to the next-to- leading order result. 
The histograms can also be turned into dynamical K-factors, which we display in the 
lower panels. The small size of the corrections to the total cross section is reflected 
only in the angular distributions, where we can see positive corrections of the order 
of 5% — 10%. Both distributions of m t j and px t get sizeable negative corrections for 
large values of these observables. For the m t i distribution, corrections reach —30% 
which has to be compared with positive +25% corrections close to the tt threshold. 
The pT t distribution is corrected down to —40% at the tails and +20% for small 
values of pT t - Overall, this leads to a distortion of the differential distributions up 
to 55% — 60%. Given that top-quark pair production at high scale is an ideal tool 
to search for various models of physics beyond the Standard Model with new gauge 
bosons like e.g. Z' , it is clear that a precise knowledge of the higher order corrections 
in this region is of significant importance. 

In Figure 0, the 6-jet kinematics is presented, where differential cross section 
distributions as a function of the averaged transverse momentum, pT b , and averaged 
rapidity, y b , of the b- and anti-6-jet are presented together with the AR bb separation. 
Both angular distributions, y b and AR bb , exhibit small positive corrections 5% — 10%, 
however, for the px b distribution we observe large and positive corrections of the order 
of +30% at the begin of the spectrum and negative of the order of —20% around 200 
GeV. 

A similar situation is observed for the W boson kinematics which is shown in 
Figure §, where the differential cross section distributions as function of the averaged 
transverse momentum px w of the W ± bosons together with an averaged rapidity yw 
of the W ± bosons are depicted. Yet again, small positive corrections of 5% — 10% 
are acquired for angular distributions as well as for low values of px w , while the tail 
of the pt w differential distribution exhibits negative corrections down to —30%. 

Subsequently, in Figure |^, differential cross section distributions as function of 
the averaged transverse momentum pT t and averaged rapidity yi of the charged lep- 
tons together with PT misa and the separation Ai?« are shown. Also here, a distortion 
of the pt ( differential distribution up to 40% is reached, while for PT misa up to 15%. 
For the angular distributions, moderate corrections up to +10% are obtained. 

And finally, in Figure the differential cross section distribution as function of 



- 19 - 



the total transverse energy defined as 



Ht = Pn + Pr- b + Pr e+ + Pt^ + PT miaa (3.8) 

is presented. In this case we observe a distortion of the differential distribution up 
to 70% - 80%. 

Overall, we can say that at the TeVatron, employing a fixed scale ft = m t , the 
NLO corrections to transverse momentum distributions are moderate. However, they 
do not simply rescale the LO shapes, but induce distortions at the level of 15% — 80%, 
which redistribute events from larger to smaller transverse momenta. The same 
applies to the invariant mass distribution of the tt pair. As for angular distributions 
we observe positive and rather modest corrections of the order of 5% — 10%. 



3.3 Results for the LHC 



Table H shows the integrated cross sections at the LHC with yfs = 7 TeV, for two 
choices of the a max parameter and for three different jet algorithms. At the central 
scale value, the full cross section receives NLO QCD corrections of the order of 



47%. Figure |TT| presents the dependence of the integrated LO cross section on the 
renormalization and factorization scales where fi = fiR = tip = £m t . The variation 
range from fi = m t /8 to ft = 8m t . In contrast to the TeVatron, the gg channel 
comprises about 76% of the LO pp cross section, followed by the qq channel with 
about 24%. In the right panel of Figure [II], the scale dependence of the NLO cross 
section is shown together with the LO one. Comparing the LO and NLO predictions, 
we find again that the large scale dependence of about +37% and —25% in the LO 
cross section is considerably reduced, down to +4% and —9% when varying the scale 
down and up by a factor 2, after including the NLO corrections. 



Algorithm 


(Xlo [fb] 


<lo =1 [fb] 


<m =om [fb] 


anti-kr 


550.54 ± 0.18 


808.46 ± 0.98 


808.29 ± 1.04 




550.54 ± 0.18 


808.67 ± 0.97 


808.86 ± 1.03 


C/A 


550.54 ± 0.18 


808.74 ± 0.97 


808.28 ± 1.03 



Table 2: Integrated cross section at LO and NLO for pp — > e Jr v e fi~v^bb + X production 
at the LHC with yfs = 7 TeV, for three different jet algorithms, the anti-kT, kx and for 
the Cambridge/Aachen jet algorithm. The two NLO results refer to different values of the 
dipole phase space cutoff a max . The scale choice is fiR = iif = mt- 
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Figure 11: Scale dependence of the LO cross section with the individual contributions of 
the partonic channels (left panel) and scale dependence of the LO and NLO cross sections 
(right panel) for the pp — > e + u e ^~h , iJ bb + X process at the LHC with ^fs = 7 TeV, where 
renormalization and factorization scales are set to the common value /j, = [ir = \xf = £mt. 



Algorithm 


a LO [fb] 


^nEo =1 [fb] 


^nEo = ° 01 [fb] 


anti-kT 


1394.72 ± 0.75 


1993.3 ± 2.5 


1993.9 ± 2.7 




1394.72 ± 0.75 


1995.2 ± 2.5 


1994.3 ± 2.7 


C/A 


1394.72 ± 0.75 


1995.0 ± 2.5 


1994.3 ± 2.7 



Table 3: Integrated cross section at LO and NLO for pp —> e + v e ^~v^bb + X production 
at the LHC with yfs = 10 TeV ; for three different jet algorithms, the anti-kT, and for 
the Cambridge /Aachen jet algorithm. The two NLO results refer to different values of the 
dipole phase space cutoff ' a ma x- The scale choice is fiR = jif = mt- 



In order to quantify the size of the non-factorizable corrections for the LHC, we 
analyze once more the narrow-width limit of our calculation following the procedure 
described in Section |3]2|. Our results are presented in Figure 12), where the depen- 
dence of the total NLO cross section together with its individual contributions, real 
emission part and the LO plus virtual corrections, are shown. Also in this case, the 
behavior is compatible with a logarithmic dependence on T t which cancels between 
real and virtual corrections. Going from NWA to the full result changes the cross 
section no more than —1.2% for our inclusive setup, which is within the expected 
uncertainty of 0(T t /m t ) of the NWA approach. 
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Figure 12: Dependence of the NLO cross section, ctt, (red solid line) and the individual 
contributions, the real emission part, cxr, (green dashed line) and the LO plus virtual part, 
cry, (blue dotted line), on the rescaling parameter C, defined as ^rescaled 

= (T t for the 

pp — > e + v e ^L~D^bb + X process at the LHC with yfs = 7 TeV. Act is defined as follows: 
Aa,(C) = (ffi(C) - <7»(C = 1))MKC = 1) with i = V,R,T. 

In Table [J[ the integrated cross sections at the LHC with \/s = 10 TeV are 
presented, once more for two choices of the a max parameter and for the three different 
jet algorithms. In this case, at the central scale value, the full cross section receives 
NLO QCD corrections of the order of 43%. 

In a next step, we compare our NLO integrated cross section with the value 
cnlo = 2097 fb for y/s = 10 TeV presented in Ref. [8]. We observe a 5% discrep- 
ancy which can perfectly be explained using the same arguments as in the TeVatron 
case, namely the on-shell top and W boson approximation applied in [8] and small 
differences between individual setups. 

As in the case of TeVatron, we have also made a comparison with Mcfm. We 
have obtained cx LO = (563.01 ± 0.63) fb, o~ nlo = (838.98 ± 1.68) fb for ^/l = 7 
TeV and a LO = (1421.05 ± 1.59) fb, a NLO = (2046.9 ± 4.3) fb for = 10 TeV, 
which constitutes a difference of 2% at LO for both cases and a difference of 4% and 
3% respectively at NLO. Moreover, both NLO results remain within our theoretical 
uncertainty of 9%, which is due to scale variation. 

Top quark production at the LHC is forward-backward symmetric in the labo- 
ratory frame as a consequence of the symmetric colliding proton-proton initial state. 
Therefore, we turn our attention to the size of NLO QCD corrections to the differ- 
ential distributions at the LHC. 

We present the differential distributions only for the y/s = 7 TeV case. In 
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Figure 13: Differential cross section distributions as a function of the invariant mass m t j 
of the top-anti-top pair, rapidity y t f of the top-anti-top pair, averaged transverse momentum 
PT t of the top and anti-top and averaged rapidity yt of the top and anti-top for the pp —> 
e + v e ^i~v^hb + X process at the LHC with yfs = 7 TeV. The blue dashed curve corresponds 
to the leading order, whereas the red solid one to the next-to-leading order result. The lower 
panels display the differential K factor. 

Figure |H| differential distributions of the tt invariant mass, m t i, together with the 
rapidity distribution, y t t, of the top-anti-top system as well as the averaged transverse 
momentum, p Tt , and the averaged rapidity y t of the top and anti-top are depicted. 
Distributions become harder in and in the invariant mass of the tt pair, moving 
from the TeVatron to the LHC expected from the higher scattering energy. 

Rapidity distributions of the tt pair and the t quark, on the other hand, get broadened 
in this transition. NLO QCD corrections to these differential distributions are always 
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Figure 14: Differential cross section distributions as a function of the averaged transverse 
momentum px b of the b-jet, averaged rapidity y b of the b-jet and Ai?^ separation for the 
pp —7- e^Veii^v^bb + X process at the LHC with yfs = 7 TeV. The blue dashed curve 
corresponds to the leading order, whereas the red solid one to the next-to-leading order 
result. The lower panels display the differential K factor. 



positive and below 50% — 60%. In case of rapidity distributions this applies for events 
concentrated within \y t i\ < 2 and \y t \ < 2 regions. 

In Figure [14], the 6-jet kinematics is presented again, but this time in the frame- 
work of the LHC. In particular, differential cross section distributions as function 
of the averaged transverse momentum, pT b , and averaged rapidity, y^, of the b- and 
anti-6-jet are presented together with the AR b i separation. Also in this case, the 
PT b distribution is harder than at the TeVatron and the y& distribution is broader. 
Clearly, the distributions show the same large and positive corrections, which turn 
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Figure 15: Differential cross section distributions as a function of the averaged transverse 
momentum pt w of the bosons and averaged rapidity yw of the bosons for the 
pp — > e + u e fi~i' fl bb + X process at the LHC with ^fs = 7 TeV. The blue dashed curve 
corresponds to the leading order, whereas the red solid one to the next-to-leading order 
result. The lower panels display the differential K factor. 

out to be relatively constant. Only in case of AR b i, corrections lead to a distortion 
of the differential distributions up to 30%. 

The W ± boson kinematics is shown in Figure [15[ where the differential cross 
section distributions as a function of the averaged transverse momentum px w of the 
W ± bosons together with an averaged rapidity yw of the W ± bosons are depicted. 
Large positive corrections of 50% — 60% are acquired for pt w differential distribution 
and rapidity distribution with events concentrated within \yw±\ < 2. The tails of 
the y w ± distribution acquire even higher NLO QCD corrections. 

Subsequently, in Figure [TBI, differential cross section distributions as function 



of the averaged transverse momentum px t and averaged rapidity yi of the charged 
leptons together with pr miss and ARu separation are shown. A small distortion of 
the px e differential distribution up to 25% is reached, while for pr misa a distortion 
up to 70% — 80% is visible. For the ye distribution, large positive and rather con- 
stant corrections up to 50% are obtained, and for the tails of the ARu distribution 
corrections of 80% — 90% are obtained. 



Finally, in Figure [17] the differential cross section distribution as function of 
the total transverse energy defined in ( p.8|) is presented. In this case we observe a 
distortion of the differential distribution up to 40%. 

Generally, we can say that for a fixed scale /i = nit at LHC, the NLO QCD 
corrections are always positive and large, at the level of 50% — 60%. Furthermore, 
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Figure 16: Differential cross section distributions as a function of the averaged transverse 
momentum pt ( of the charged leptons, averaged rapidity yi of the charged leptons, PT misa 
and ARi£ for the pp — > e + v e ii~v^bb + X process at the LHC with ^/s = 7 TeV. The blue 
dashed curve corresponds to the leading order, whereas the red solid one to the next-to- 
leading order result. The lower panels display the differential K factor. 

they are relatively constant. Exceptions are the rapidity distributions, which are 
only constant in the central region, and the PT miss and Ht distributions, which are 
distorted up to 40% - 80%. 



4. Conclusions 

In this paper, we have presented, for the first time, a computation of the NLO QCD 
corrections to the full decay chain pp(pp) — > tt — > W + W~bb — > e + v e ^~v fJ bb + X. 
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Figure 17: Differential cross section distribution as a function of the total transverse 
energy, Ht, for the pp — > e + v e ^~v fJ bb + X process at the LHC with yfs = 7 TeV. The 
blue dashed curve corresponds to the leading order, whereas the red solid one to the next- 
to-leading order result. The lower panels display the differential K factor. 

All off-shell effects of top quarks and W gauge bosons have been included in a fully 
differential manner which allows us to compute an arbitrary observable in terms of 
jets, charged leptons and missing transverse energy within experimentally relevant 
selection criteria with NLO QCD accuracy In order to illustrate the capabilities 
of the program, the total cross section and its scale dependence, as well as several 
differential distributions at the TeVatron run II and the LHC have been given. More- 
over, in case of the TeVatron the forward-backward asymmetry of the top has been 
recalculated. We have found that with inclusive selection cuts, the forward-backward 
asymmetry amounts to A t FB = 0.051 ± 0.0013. Furthermore, the impact of the NLO 
QCD corrections on integrated cross sections at the TeVatron is small, of the order 
2.3%. At the LHC we have obtained NLO QCD corrections at the level of 47% and 
43% for -y/i = 7 TeV and y/s = 10 TeV respectively. A study of the scale depen- 
dence of our NLO predictions indicates that the residual theoretical uncertainty due 
to higher order corrections is 8% for the TeVatron and 9% for the LHC. 
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Note added 

Independently of our calculation, another group has evaluated NLO QCD corrections 
to WWbb production with leptonic decays of gauge bosons, and has presented them 
in [83]. We have cross checked the results, applying the narrow width approximation 
for the W bosons as in that publication, and have obtained perfect agreement for 
integrated LO and NLO cross sections within statistical errors. 
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